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California,  Los  Angeles,  California  90089,  Department  of  Chemistry,  University  of  Warwick, 
Coventry  CV4  7AL,  West  Midlands,  UK,  and  Department  of  Chemistry,  Iowa  State  University, 
Ames,  Iowa 

Abstract 

Ab  initio  molecular  orbital  theory  has  been  used  to  calculate  accurate  enthalpies  of  formation  and 
adiabatic  electron  affinities  or  ionization  potentials  for  N3,  N3 ,  Ns^,  and  N5  from  total 
atomization  energies.  Bom-Haber  cycle  calculations,  using  estimated  lattice  energies  and  the 
adiabatic  ionization  potentials  of  the  anions  and  electron  affinities  of  the  cations,  permit  for  the 
first  time  reliable  stability  predictions  for  the  hypothetical  N5^3  and  Ns^s  salts.  Contrary  to 
previous  predictions,  it  is  shown  that  neither  salt  can  be  stabilized,  decomposing  spontaneously 
to  N3  radicals  and  N2.  This  conclusion  was  experimentally  confirmed  by  low-temperature 
metathetical  reactions  between  NsSbFe  and  alk^i  metal  azides  in  different  solvents,  resulting  in 
violent  reactions  with  spontaneous  nitrogen  evolution.  It  is  also  shown  that  the  vertical  ionization 
potentials  and  electron  affinities,  used  in  previous  stability  predictions,  must  not  be  used  in 
stability  predictions  involving  species,  such  as  N5  ions,  because  the  resulting  N5  radicals  are  not 
vibrationally  stable  and  undergo  spontaneous  further  decomposition.  The  use  of  the  vertical 
values,  that  neglects  this  highly  exothermic  secondary  decomposition  of  N5  to  N3  and  N2,  results 
in  large  errors  of  about  100  kcal/mol  for  each  N5  unit.  It  is  also  shown  that  a  previous 
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comparison  of  the  density  and  energy  density  of  hydrazine  with  those  estimated  for  Ns'^s'  is 
significantly  in  error. 
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Introduction 

The  discovery  of  stable  Ns'^  salts, the  long-known  existence  of  stable  Ns'  salts, ^  and  the 
recent  experimental  detection  of  the  N5'  anion'*  have  stimulated  research  in  polynitrogen 
chemistry  and  a  search  for  nitrogen  allotropes.  The  high  endothermicity  of  polynitrogen 
compounds  renders  their  preparation  and  handling  very  difficult,  and  reliable  stability  predictions 
are  important  for  the  selection  of  suitable  target  compounds.  Unfortunately,  only  the  enthalpy  of 
formation  of  N2  is  well  established.  The  enthalpies  of  formation  of  other  crucial  species,  such  as 
N3,  Ns',  N5,  Ns"^  and  Ns',  are  not  accurately  known.  There  have  been  a  number  of  theoretical 
studies  on  allotropes  of  nitrogen  using  ab  initio  molecular  orbital  theory  without  the  use  of 
empirical  parameters,  but  these  studies  have  not  reliably  established  the  energetics  of  these 
compounds.  In  particular  we  note  the  important  work  of  the  Bartlett  group  on  these  species,^’®  the 
work  of  Nguyen  and  Ha  on  decomposition  mechanisms  for  N5’‘,  x=  -1, 0,  +1,^  and  that  of  Martin 
and  co-workers  on  N3.® 

We  have  developed  an  approach^  to  the  reliable  calculation  of  molecular  thermodynamic 
properties,  notably  enthalpies  of  formation,  based  on  ab  initio  molecular  orbital  theory  without 
reliance  on  embedded  empirical  parameters.  Our  approach  involves  an  accurate  calculation  of 
the  total  atomization  energy  of  a  molecule  and  combines  this  value  with  known  enthalpies  of 
formation  of  the  atoms  to  calculate  the  molecular  enthalpy  of  formation  at  0  K.  This  method 
starts  with  coupled  cluster  theory,  including  a  perturbative  triples  correction  (CCSD(T)),*°’”’*^ 
combined  with  the  correlation-consistent  basis  sets*^’*'*  extrapolated  to  the  complete  basis  set 
(CBS)  limit  to  treat  the  correlation  energy  of  the  valence  electrons.  This  step  is  followed  by  a 
number  of  smaller  corrections  that  are  presumed  to  be  additive,  including  core-valence 
interactions  and  relativistic  effects,  both  scalar  and  spin-orbit.  Finally,  one  must  include  the  zero 
point  energy  obtained  either  from  experiment,  theory,  or  some  combination.  The  standard 
enthalpies  of  formation  of  compounds  at  298  K  can  then  be  calculated  from  the  0  K  value  and 
other  calculated  properties  by  using  standard  thermodynamic  and  statistical  mechanics 
expressions.*^ 

Another  important  problem  involves  the  reliable  prediction  of  whether  the  combination 
of  a  stable  pol3mitrogen  cation,  such  as  Ns"^,  with  a  stable  polynitrogen  anion,  such  as  N3*,  can 
lead  to  a  stable  ionic  nitrogen  allotrope,  such  as  Ns'T^s'.  Whereas  such  predictions  can  be  made 
with  the  help  of  suitable  Bom-Haber  cycles,  they  require  the  availability  of  reliable  values  for  the 
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lattice  energies  of  the  salts,  the  first  ionization  potential  of  the  anion  and  the  electron  affinity  of 
the  cation,  and  an  excellent  knowledge  of  the  nature  and  stability  of  the  most  likely 
decomposition  products.  Although  several  theoretical  studies  on  the  stability  of  Ns’^s  and 
Ns'l^s'  have  been  published,^’®’^^  none  of  these  have  used  appropriate  Bom-Haber  cycles  and 
input  data.  Therefore,  the  most  important  question  has  not  been  answered  conclusively.  Is  the 
lattice  energies  of  these  salts  sufficient  to  stabilize  them  as  solids  and  prevent  a  spontaneous 
electron  transfer  from  the  anion  to  the  cation  causing  spontaneous  decomposition?  A  further 
important  question  deals  with  the  problem  as  to  whether  calculated  vertical  electron  affinities 
and  ionization  potentials  can  be  used  for  these  predictions.  For  example,  the  vertical  electron 
affinity  calculated®  for  differs  significantly  from  the  experimentally  determined  electron 
affinity,^  and  there  appears  to  be  confusion  with  regard  to  which  value  should  be  used  for  the 
stability  predictions.  This  paper  provides  conclusive  answers  to  these  questions  and,  for  Ns'^s'^ 
provides  experimental  data  in  support  of  our  theoretical  predictions. 

Experimental  Section 

Caution!  Reactions  of  N 5^  salts  with  azide  ions  are  violent  and  can  result  in  explosions. 
Therefore,  these  materials  should  be  handled  only  on  a  small  scale  with  appropriate  safety 
precautions  (face  shield,  leather  gloves,  and  protective  clothing). 

Materials  and  Apparatus.  All  reactions  were  carried  out  in  a  Teflon-PFA  double  U-tube 
apparatus  that  consisted  of  a  reaction  U-tube,  a  porous  Teflon  filter  assembly,  and  a  receiver  U- 
tube.  Volatile  materials  were  handled  on  a  stainless  steel/Teflon-FEP  vacuum  line.^^  Nonvola¬ 
tile  solids  were  handled  in  the  dry  nitrogen  atmosphere  of  a  glove  box.  Infrared  spectra  were 
recorded  on  a  Mattson  Galaxy  FT-IR  spectrometer  using  dry  powders  pressed  between  AgCl 
windows  in  an  Econo  press  (Barnes  Engineering  Co.).  Raman  spectra  were  recorded  on  a  Bruker 
Equinox  55  FT-RA  spectrometer  using  a  Nd-YAG  laser  at  1064  nm  and  Pyrex  melting  point 
capillaries  as  sample  containers. 

The  NsSbFe  was  prepared  based  on  our  previously  reported  procedure.  The  CsNs 
(Aldrich  Chemical  Company)  was  pretreated  with  an  excess  of  (CH3)3SiN3  to  remove  any  traces 
of  moisture.  The  SO2  (Air  Products)  was  used  as  received. 

Reactions  of  NsSbFg  with  CsNs.  In  a  typical  experiment,  N5SbF6  (2.285  mmol)  and  CSN3  (2.267 
mmol)  were  loaded  in  the  dry  box  into  a  Teflon-PFA  reaction  vessel  that  was  closed  by  a  Teflon 
valve.  On  the  vacuum  line,  anhydrous  SO2  (5.0  ml,  liquid)  was  added  at  -196  °C.  The  reaction 


6 


mixture  was  warmed  to  the  melting  point  of  SO2  (-64  °C).  After  the  SO2  had  partially  melted  and 
the  slurry  of  reactants  was  gently  agitated,  a  violent  reaction  took  place  producing  a  bright 
yellow  flame.  The  reactor  was  quickly  quenched  with  liquid  nitrogen.  A  check  for 
noncondensible  gases  at  -196  °C  revealed  9.05  mmol  of  N2  as  expected  for  the  complete 
decomposition  of  2.267  mmol  N5N3  to  N2.  Pumping  off  all  volatiles  at  -64  °C  produced  0.8890  g 
of  a  white  solid  residue  (weight  expected  for  2.267  mmol  CsSbFe  =  0.836  g).  The  extra  weight 
was  due  to  some  residual  SO2.  Raman  and  infrared  spectra  of  the  white  solid,  after  pumping  at 
room  temperature,  showed  it  to  be  pure  CsSbFg. 

This  reaction  can  be  moderated  by  the  separate  dissolution  of  the  starting  materials  in 
SO2  and  the  combination  of  the  resulting  clear  solutions  at  -64  °C.  Under  these  conditions,  the 
decomposition  reaction  proceeded  rapidly  but  without  violence  and  gave  the  same  results  as 
described  above. 

Reactions  of  NsSbFg  with  CsNa.  To  equimolar  amounts  of  NaNs  and  NsSbFe,  CHF3  was  added 
as  a  solvent  at  -196  °C,  and  the  resulting  mixture  was  warmed  to  -95  and  then  to  -64  °C.  No 
N2  evolution  was  observed  under  these  conditions.  The  CHF3  was  pumped  off  at  -95  °C  and  SO2 
was  added  as  a  solvent  at  -196  °C.  The  mixture  was  warmed  to  -64  °C,  but  again  no  N2 
evolution  was  observed.  Pumping  off  the  SO2  at  -64  “C  resulted  in  a  white  solid  residue  that, 
based  on  its  room  temperature  Raman  spectrum,  consisted  of  a  mixture  of  unreacted  NaNs  and 
NsSbFe,  indicating  that  no  reactions  had  occurred  due  to  the  very  low  solubility  of  NaN3  in  these 
solvents. 

Computational  Approaches 

Enthalpies  of  Formation.  The  augmented  correlation  consistent  basis  set  aug-cc-pVnZ  (n  =  D, 
T,  was  used  for  nitrogen.  Only  the  spherical  components  (5-J,  l-f  9-g  and  1  l-h)  of  the 
Cartesian  basis  functions  were  used.  This  family  of  basis  sets  contains  an  additional  shell  (e.g. 
spd  for  DZ)  of  diffuse  functions  that  are  necessary  for  an  accurate  description  of  anions. 
Calculations  were  performed  using  the  MOLPRO,^®  Gaussian,^®  and  GAMESS^’  suites  of 
programs.  The  open-shell  CCSD(T)  calculations  were  ckried  out  at  the  R/UCCSD(T)  level.  In 
this  approach,  a  restricted  open  shell  Hartree-Fock  (ROHF)  calculation  was  initially  performed 
in  order  to  generate  the  set  of  molecular  orbitals  and  the  spin  constraint  was  relaxed  in  the 
coupled  cluster  calculation.^^’^^’^'^  The  CCSD(T)  total  energies  were  extrapolated  to  the  CBS 
limit  by  using  a  mixed  exponential/Gaussian  function  of  the  form: 
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E(n)  =  EcbS  +  A  exp[-(n-l)]  +  B  exp[-(n-l)2]  (1) 

with  n  =  2  (aug-cc-pVDZ),  3  (aug-cc-pVTZ),  etc.,  as  first  proposed  by  Peterson  et  al.^^  This 
extrapolation  method  has  been  shown  to  yield  atomization  energies  in  the  closest  agreement  with 
experiment  by  a  small  measure  as  compared  to  other  extrapolation  approaches  up  through  n  = 
4.* 

Most  geometries  were  optimized  at  the  frozen  core  CCSD(T)  level  with  various 
augmented  correlation-consistent  basis  sets.  With  the  exception  of  the  Ns  radical,  all  geometries 
were  optimized  with  the  aug-cc-pVTZ  or  aug-cc-pVQZ  basis  sets.  If  the  CCSD(T)/aug-cc-pVQZ 
geometry  was  not  available,  the  CCSD(T)/aug-cc-pVTZ  geometry  was  used  for  the 
CCSD(T)/aug-cc-pVQZ  calculation.  Geometries  for  the  Ns  radicals  were  obtained  from 
ZAPT2/6-311+G(2df)  and  CCSD(T)/aug-cc-pVpZ  calculations.  Vibrational  frequencies  were 
calculated  at  the  CCSD(T)  and  MP2  levels  with  the  aug-cc-pVDZ  basis  set.  Calculations  were 
performed  on  SGI,  IBM  and  Apple  G4  computers. 

Core-valence  corrections,  AEcv,  were  obtained  at  the  CCSD(T)/cc-pCVTZ  level  of 
theory. Scalar  relativistic  corrections  (AEsr),  which  account  for  changes  in  the  relativistic 
contributions  to  the  total  energies  of  the  molecule  and  the  constituent  atoms,  were  included  at  the 
CI-SD  (configuration  interaction  singles  and  doubles)  level  of  theory  using  the  cc-pVTZ  basis 
set.  AEsr  is  taken  as  the  sum  of  the  mass-velocity  and  1-electron  Darwin  (MVD)  terms  in  the 
Breit-Pauli  Hamiltonian.^^  Because  N  has  a  ^S  ground  state,  there  are  no  atomic  spin  orbit 
corrections  to  the  total  atomization  energies.  The  spin  orbit  coupling  constant  for  N3  is  71.3  cm'^ 
based  on  experiment^*  leading  to  a  spin  orbit  correction  of  0.10  kcal/mol. 

By  combining  our  computed  XPo  values  with  the  known  enthalpies  of  formation^  at  0  K 
for  N  of  AHf°(N)  =  1 12.53  +  0.02  kcal  mor\  we  can  derive  AH/  values  for  the  molecules  under 
study  in  the  gas  phase.  We  obtain  enthalpies  of  formation  at  298  K  by  following  the  procedures 
outlined  by  Curtiss  and  coworkers^^  and  by  using  the  stationary  electron  convention,  i.e.,  we  do 
not  assume  a  change  in  the  enthalpy  of  formation  of  the  electron  as  a  function  of  temperature. 
Lattice  Energies.  For  this  study,  we  employ  the  equation^®: 

UL  =  2/[aV;„-''*  +  P]  (2) 

to  estimate  the  lattice  energy,  Ul,  of  the  salts  Ns'^s'  and  N5''N3'  from  estimated  ion  volumes.  I  is 
the  ionic  strength  (=  1),  Vm  is  the  molecular  (formula  unit)  volume  of  the  lattices  involved  which 
is  equal  to  the  sum  of  the  individual  ion  volumes  of  the  cation,  V+  and  anion,  V.  and  a  and  P  take 
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the  values  28.0  kcal  moF^  nm  and  12.4  kcal  moF’  respectively  for  1:1  salts.  The  individual  ion 
volumes  can  be  estimated  from  an  ion  volume  database,  inferred  in  some  cases  from 
established  crystal  structure  data,  or  calculated.  In  order  to  calculate  the  volumes  of  the  ions,  we 

o  1 

chose  to  use  the  volumes  that  we  have  used  in  free  energy  of  solvation  calculations.  The 
electron  densities  were  calculated  at  the  B3LYP/6-31+G*  level^^  and  the  volume  was  taken  to  be 
that  inside  the  0.001  a.u.  contour  of  the  electron  density. 

Results  and  Discussion 
Calculated  Geometries 

The  calculated  geometries  for  the  molecules,  Na,  N3,  Ns",  N5',  and  Ns"^  are  summarized  in 
Table  1. 

Nj:  As  expected,  the  geometry  for  N2  at  the  CCSD(T)/aug-cc-pVQZ  level  is  in  excellent 
agreement  with  the  experimental  value.^^  Theory  overestimates  experiment  by  0.0028  A.  For 
comparison  purposes,  we  also  optimized  the  N-N  bond  length  with  the  larger  aug-cc-pV5Z  basis 
set.  This  produced  a  small  0.0010  A  decrease  in  the  bond  length,  indicating  that  aug-cc-pVQZ 
bond  lengths  should  be  within  approximately  0.001  A  of  the  basis  set  limit  for  the  other 
molecules  that  could  be  treated  at  this  level  of  theory.  If  the  frozen  core  approximation  is 
removed  via  calculations  with  the  core/valence  (CV)  cc-pCVQZ  basis  set,  another  0.0020  A 
shortening  of  the  bond  length  is  observed.  The  final  CCSD(T)/aug-cc-pV5Z  +  CV  bond  length 
of  1.0975  A  is  almost  identical  to  the  experimental  value  at  1.0977  A 

N3:  The  geometry  for  N3  (D^h)  is  in  excellent  agreement  with  the  experimental  value.^ 
The  calculations  predict  a  value  for  the  NN  bond  that  is  slightly  shorter  than  the  experimental 
value  in  contrast  to  N2  where  the  calculations  predict  a  longer  bond  at  the  CCSD(T)/aug-cc- 
pVQZ  level. 

Na’:  The  calculated  geometry  for  Ns'  is  also  in  very  good  agreement  with  experiment, 
with  the  calculated  NN  bond  length  again  being  slightly  shorter  than  experiment.^^ 

Ns"^:  The  geometry  for  Ns"^  has  been  measured  in  the  crystal  by  x-ray  diffraction 
techniques.  The  calculated  values  are  in  good  agreement  with  the  crystal  structure  values 
considering  the  uncertainty  in  the  experimental  values  and  the  different  environment.  The 
calculated  geometrical  parameters  for  Ns"^  are  the  same  as  those  reported  by  Fau  and  Bartlett  at 
the  CCSD(T)/aug-cc-pVTZ  level,^  whereas  the  geometries  calculated  for  Ns'  at  the  same  level 
differ  slightly. 
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Ns':  The  value  of  1.334  A  calculated  for  the  NN  bond  length  in  Ns'  at  the  CCSD(T)/6- 
311+(G(d,p)  level^  lies  between  our  values  calculated  at  the  CCSD(T)  level  with  the  aug-cc- 
pVDZ  and  aug-cc-pVTZ  basis  sets,  respectively. 

Ns  (open  isomer):  ZAPT2^^  calculations  with  the  6-311+G(2df)  basis  set  followed  by 
CCSD(T)/cc-pVrZ  single  point  calculations  confirm  previous  conclusions^  that  the  Ns  radical 
starting  from  this  geometry  is  vibrationally  unstable.  Starting  from  the  optimized  Ns^  structure, 
the  neutral  species  dissociates  without  a  barrier  to  N3  +  N2  as  shown  in  Figure  1. 

Ns  (cyclic  isomer):  The  Ns  radical  cyclic  isomer  would  be  generated  by  detaching  an 
electron  from  cyclic  Ns'.  The  HOMO  of  Ns'  is  a  degenerate  in-plane  orbital  involving  essentially 
lone  pair  interactions.  Production  of  the  Ns  radical  leads  to  (e)^  occupancy  which  Jahn-Teller 
distorts  to  a  Civ  structure.  Two  states  can  be  derived  from  the  (e)  occupancy,  a  Ai  state  and  a 
state(a  state  would  be  a  n  radical).  The  geometry  for  the  cyclic  Ns  (Civ)  ^Ai  radical  was 
found  to  be  a  stationary  point  on  the  ZAPT2/6-311+G(2df)  potential  energy  surface.  The  cyclic 
Ns  radical  at  the  ZAPT2/6-311+G(2df)  geometry  is  66.8  kcal/mol  above  the  N3  +  Ni  asymptote 
at  the  CCSD(T)/cc-pVTZ  level.  The  longest  NN  ZAPT/6-311+G(2df)  bond  distance  is  that  in 
the  base  of  the  pentagon,  1.354  A.  The  shortest  NN  bonds  are  those  connected  to  the  apex,  1.292 
A.  We  took  this  geometry  and  initiated  an  optimization  at  the  R/UCCSD(T)/aug-cc-pVDZ  level 
for  both  the  ^Ai  and  ^Bi  states  under  the  constraint  of  Civ  symmetry.  The  ^Bi  state  smoothly 
dissociated  to  the  N3  +  Ni  asymptote  showing  that  on  the  CCSD(T)/aug-cc-pVDZ  surface,  this 
state  is  not  a  minimum.  The  N3  distance  converges  faster  than  the  Ni  distance.  The  Ai  state  is 
difficult  to  optimize  but  finally  optimizes  to  a  structure  with  2  elongated  Ni  molecules 
complexed  weakly  to  an  N  atom  (see  Table  1)  as  it  cannot  go  to  the  products  Ni  +  N3.  At  the 
CCSD(T)/aug-cc-pVTZ  level  at  the  optimized  CCSD(T)/aug-cc-pVDZ  geometry,  the  ^Ai  state  is 
60.2  kcal/mol  above  the  Ni  +  N3  asymptote.  These  calculations  strongly  suggest  that  on  electron 
detachment  from  N5',  the  resulting  radical  has  a  path  on  the  ^Bi  surface  that  leads  directly  to  the 
N3  -I-  Ni  asymptote  and  that  there  is  no  stable  cyclic  N5  radical  structure  as  already  noted 
previously.^  Frenking  and  co-workers  as  part  of  a  study  of  FeCri^-Ns)!  report  a  structure  for  the 
^A]  state  of  the  N5  radical  at  the  B3LYP/6-31G(d)  level  which  could  be  a  minimum  although  the 
structure  is  not  discussed. 

Calculated  Frequencies 
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The  calculated  vibrational  frequencies  for  N2,  N3,  Ns’,  Ns’,  Ns  (cyclic)  and  Ns'^  are 
summarized  in  Table  2. 

N2:  The  harmonic  frequency  at  the  frozen  core  CCSD(T)/aug-cc-pVQZ  level  is  within  5 
cm*’  of  the  experimental  value/^  Even  at  the  CCSD(T)/aug-cc-pVTZ  level,  the  N2  frequency  is 
within  20  cm’’  and  is  within  40  cm’’  at  the  CCSD(T)/aug-cc-pVDZ  level.  This  level  of 
agreement  is  in  part  fortuitous.  Whereas  the  aug-cc-pV5Z  basis  set  frequency  is  identical  to  the 
experimental  value,  introduction  of  core/valence  correlation  increases  (Og  by  9.2  cm’’.  Thus,  the 

final  CCSD(T)/aug-cc-pV5Z  +  CV  value  is  9  cm’’  larger  than  experiment. 

N3:  For  Ns,  the  frequencies  calculated  at  the  MP2/aug-cc-pVDZ  level  are  in  qualitative 
agreement  with  the  experimental  values^*  for  the  7t  bend  and  the  symmetric  stretch  and  in  good 
agreement  for  the  antisymmetric  stretch.  In  this  case,  we  have  used  the  experimental  frequencies 
to  calculate  the  zero  point  energy. 

Ns':  For  Ns‘,  the  calculated  antisymmetric  harmonic  stretch  is  in  excellent  agreement  with 
the  gas  phase  experimental  value.^'’’^^  The  value  for  the  V2  bending  frequency  for  this  ion  has 
been  measured  in  the  solid  and  a  free  space  estimate  for  Vj  has  been  obtained  based  on  the 
spectrum  measured  in  the  solid.^’'”  For  this  ion,  we  used  the  experimental  frequencies  for  the 
zero  point  energy  determination.  Use  of  the  theoretical  frequencies  introduces  a  difference  of 
0.19  kcal/mol  because  the  calculated  bend  and  antisymmetric  stretch  are  ~  50  cm’’  lower  than  the 
experimental  values.  The  frequencies  calculated  by  us  for  Vj  and  V3  are  in  reasonable  agreement 
with  those  previously  calculated  at  the  CASPT2  level’^  and  at  the  CEPA  level.'” 

Ns’;  For  this  ion,  we  used  the  frequencies  calculated  at  the  CCSD(T)/aug-cc-pVDZ  level 
for  the  zero  point  energy  calculation  as  there  are  no  experimental  values.  Quite  good  agreement 
is  found  for  the  frequencies  of  N5’  calculated  at  the  MP2/aug-cc-pVDZ  level  and  at  the 
CCSD(T)/aug-cc-pVDZ  level.  The  CASPT2  frequencies’®  are  similar  to  our  values. 

Ns  (cyclic  isomer):  Frequencies  were  calculated  at  the  ZAPT2/6-311+G(2df)level  of 
theory  for  the  ^Ai  state.  No  experimental  values  are  available  for  this  species.  As  noted  above, 
the  ^Bi  stmcture  is  not  a  minimum  for  the  cyclic  radical  at  the  R/UCCSD(T)/aug-cc-pVDZ  level. 

Ns"^:  For  Ns"^,  there  is  reasonable  agreement  between  the  frequencies  calculated  at  the 
MP2/aug-cc-pVDZ  level  and  the  average  of  the  frequencies  taken  from  measurements  in  the 
solid  state.^  The  largest  discrepancy  is  found  for  the  highest  frequency  bs  mode  with  a  difference 
of  -140  cm’’.  The  calculated  as  and  bi  modes  deviate  by  -50  cm’’  from  the  experimental  values. 
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The  fact  that  the  differences  cancel  out  leads  to  calculated  and  experimental  zero  point  energies 
that  are  essentially  identical.  We  used  the  average  experimental  values  in  our  zero  point  energy 
calculations.  Our  calculated  frequencies  are  in  qualitative  agreement  with  those  calculated  at  the 
CCSD(T)/6-311-fG(2d)  and  B3LYP/6-311+G(2d)  levels.*’^  The  CCSD(T)  values  are  in  general 
lower  than  our  MP2  values  except  for  the  a2  and  b]  modes  which  are  higher.  The  CASPT2 
frequencies  differ  from  ours  more  than  might  be  expected,  especially  for  the  two  highest 
modes.’® 

Calculated  Enthalpies  of  Formation 

The  total  atomization  energies  are  given  in  Table  3.  The  AH/*  values  were  converted  to 
298  K  following  the  procedure  of  Curtiss  and  coworkers.’® 

Na:  As  shown  in  Table  3,  the  error  in  the  atomization  energy  of  N2  is  0.76  kcal  mol’’ 
when  compared  to  experiment;^^  This  leads  to  an  error  in  the  calculated  enthalpy  of  formation  at 
0  K  of  0.76  kcal  mol'’.  We  can  thus  estimate  that  the  error  bars  in  our  heats  of  formation  are 
likely  to  be  on  the  order  of  +  1  kcal/mol. 

N3:  The  enthalpy  of  formation  of  N3  has  been  reported  by  a  number  of  authors.  The 
NIST  value^^  of  100  +  5  kcal  mol'’  at  0  K  is  almost  10  kcal  mol'’  below  our  calculated  value. 
Our  calculated  value  of  109.3  kcal  mol  ’  is  in  much  better  agreement  with  the  value  of  112  +  5 
based  on  the  experimental  measurement  of  the  acidity  of  HN3  and  the  measurement  of  the 
electron  affinity  of  N3  Our  value  of  109.3  kcal  mol'’  is  in  excellent  agreement  with  the 
scaled  value  of  109.25  kcal  mol'’  given  by  Martin  et  al.^  The  calculated  atomization  energies 
show  that  reaction  (3) 

N3("n)  (g)  - ►  N(^S)  (g)  +  N2(’E0(g)  (3) 

is  only  4  kcal  mol'’  endothermic  at  0  K. 

N3':  The  enthalpy  of  formation  of  N3'  has  been  determined  from  gas  phase  acidity 
measurements  and  the  enthalpy  of  formation  of  HN3  to  be  48  ±  2  kcal  mol'’.'’^  This  value  is  in 
excellent  agreement  with  our  calculated  value  of  47.2  kcal  mol'’.  Fau  and  Bartlett  calculated  the 
enthalpy  of  formation  of  N3'  by  a  variety  of  methods.®  Their  best  estimated  enthalpy  of 
formation  of  47.8  kcal  mol'’  is  in  good  agreement  with  our  value,  although  we  note  that  their 
CGSD(T)/aug-cc-pVTZ  value  of  45.7  kcal  moT’  is  too  low  by  1.5  kcal  mol'’.  The  electron 
affinity  of  N3  has  been  measured  to  be  2.68  +  0.01  eV.'’®  Our  calculated  adiabatic  electron 
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affinity  of  2.69  eV  is  in  excellent  agreement  with  this  result.  With  an  ANO  (4s3p2dlf)  basis  set, 
the  CASPT2  electron  affinity  was  predicted  to  be  2.64  eV.^^ 

Ns"^:  The  enthalpy  of  formation  of  calculated  by  us,  is  353.4  kcal  mol'*.  This  result 
is  in  reasonable  agreement  with  the  value  of  347.1  kcal  mor\  calculated  by  Fau  and  Bartlett  at  a 
lower  level.^  Christe  and  co-workers*’^  predicted  the  enthalpy  of  formation  of  to  be  351.1 
kcal  mol'*  at  the  CCSD(T)  level  with  a  polarized  double  zeta  basis  set.  Nguyen  and  Ha  predict  a 
heat  of  formation  of  351.3  +  3.6  kcal/mol  based  on  CCSD(T)/6-311+  G(3df)  calculations  and  the 
experimental  energy  of  Ns.^ 

Our  value  for  Ns"^  can  be  used  to  estimate  the  adiabatic  electron  affinity  of  reaction  (4). 

N5^(g)-he--^N3(g)  +  N2(g)  (4) 

The  intermediate  Ns  radical  starting  from  the  Ns'*’  geometry  (open  structure)  is  unstable  and  will 
spontaneously  dissociate  to  the  first  stable  dissociation  products,  N3(g)  +  N2(g),  which  will  only 
form  a  weakly  bound  yan  der  Waals  complex  with  a  binding  energy  on  the  order  of  less  than  1 
kcal  mol'*.  Thus,  the  energy  of  reaction  (4)  provides  the  adiabatic  electron  affinity  of  Ns^  Our 
calculated  value  of  243.3  kcal  mol'*  (10.55  eV)  is  much  higher  than  the  vertical  electron  affinity 
of  6.04  eV  reported  by  Fau  et  al^  or  the  adiabatic  ionization  energy  of  7.4  +  0.2  eV  reported  by 
Nguyen  and  Ha,^  but  in  excellent  agreement  with  our  experimental  bracketing  study^  that 
showed  the  electron  affinity  of  Ns’*’  to  fall  between  10.52  and  11.48  eV.  In  the  calculation  of 
vertical  electron  affinities  or  ionization  potentials,  the  geometries  of  the  original  species  are 
retained  for  the  products.  Therefore,  vertical  electron  affinities  or  ionization  potentials  should 
not  be  used  for  thermochemical  calculations.  Only  when  the  direct  reaction  products  are 
vibrationally  stable  species  with  geometries  similar  to  those  of  the  starting  materials,  will  their 
vertical  EA  or  IP  values  approximate  the  adiabatic  values.  If  this  is  not  the  case,  and  the  products 
can  undergo  further  decomposition  that  does  not  involve  any  significant  activation  energy 
barrier,  then  the  decomposition  energy  to  the  first  well  defined,  vibrationally  stable  products 
must  be  added  to  the  vertical  electron  affinity  or  ionization  energy  values. 

Ns':  The  enthalpy  of  formation  of  Ns'  is  calculated  to  be  62.2  kcal  mol'*,  in  excellent 
agreement  with  the  value  of  62.1  +  3.6  kcal/mol  given  by  Nguyen  and  Ha.^  This  value  can  be 
used  to  calculate  the  adiabatic  ionization  potential  of  Ns'  in  the  same  manner  as  discussed  above 
for  Ns'*’.  As  noted  above,  loss  of  an  electron  from  the  degenerate  HOMO  of  Ns'  leads  to  either  a 
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”Bi  or  a  Ai  radical.  We  have  shown  that  formation  of  the  Bj  radical  provides  a  pathway  leading 
to  Ns  +  Na  so  the  following  channel  (Reaction  (5))  represents  the  adiabatic  ionization  process: 

Ns'  (g)  ->  Ns  (g)  +  Ns  (g)  +  e'  [47.9  kcal  mof’  (2.08  eV)]  (5) 

Use  of  this  channel  gives  an  adiabatic  ionization  potential  of  Ns'  of  2.08  eV.  This  is  much  lower 
than  the  vertical  ionization  potential  of  5.58  eV  reported  by  Fau  et.al^  or  of  5.7  eV  as  given  by 
Nguyen  and  Ha.^  Nguyen  and  Ha  do  suggest  that  the  cyclic  Ns  radical  will  spontaneously 
dissociate  to  Ns  +  Ns.  As  shown  below,  the  use  of  vertical  instead  of  the  appropriate  adiabatic 
electron  affinity  and  ionization  potential  values  in  the  Bom-Haber  cycles  dramatically  influences 
the  stability  predictions  for  species  containing  either  Ns"^  or  Ns'  ions.  For  a  salt,  such  as  Ns'*Ns’, 
containing  two  Ns  ions,  the  error  could  be  doubled  in  the  Bom-Haber  cycle  calculations 

The  important  dissociation  channel  for  Ns'  is  the  formation  of  Ns'  +  Ns.  We  predict  this 
channel  to  be  exothermic  by  14.2  kcal  mol'^  at  0  K.  This  value  is  in  good  agreement  with  the 
value  of  14.3  kcal  mol'^  reported  by  Nguyen  and  Ha  at  the  CCSD(T)/aug-cc-pVTZ  level.^  It  is 
important  to  know  if  there  is  a  barrier  separating  Ns'  from  the  asymptotic  products  Ns'  and  Ns. 
Nguyen  and  Ha  calculated  the  energy  barrier  to  dissociation  as  27.7  kcal  mol'*  at  the 
CCSD(T)/aug-cc-pVTZ  level  including  zero  point  corrections  based  on  a  geometry  calculated  at 
the  CCSD(T)/6-311-i-G(d)  level.^  We  have  calculated  the  energy  of  the  transition  state  at  the 
CCSD(T)/CBS  level  with  Nguyen  and  Ha’s  geometry  and  with  the  geometry  optimized  by  us  at 
the  MP2/aug-cc-pVDZ  level.  We  calculated  the  geometry  at  the  latter  level  in  order  to  get  a 
value  of  the  zero  point  energy  as  well  as  the  value  of  the  imaginary  frequency  characterizing  the 
transition  state.  The  results  are  given  in  Tables  1  to  3.  Our  calculated  energy  barrier  using  the 
MP2  geometry  is  29.4  kcal  mol'*  and  when  the  zero  point  energy  difference  is  included,  the 
barrier  is  27.2  kcal  mol'*.  The  energy  barrier  using  the  Nguyen  and  Ha  geometry  is  0.61  kcal 
mol'*  higher  in  energy  than'  our  value  based  on  the  MP2  geometry  showing  that  the  two 
geometries  give  essentially  identical  barrier  heights.  The  barrier  calculated’  by  Nguyen  and  Ha 
is  very  similar  to  our  CCSD(T)/CBS  value  using  their  geometry.  At  the  CASPT2/ANO(4s3p2d) 
level,  the  barrier  is  predicted  to  be  26  kcal  mol'*  and  the  dissociation  energy  to  N3'  and  N2  is 
predicted  to  be  exothermic  by  9  kcal  mol'*,*^  4  kcal  mol'*  lower  than  our  value.  A  calculation  of 
the  imaginary  frequency  for  the  transition  state  of  the  N5'  decomposition  shows  a  large  value  of 
862i  cm'*,  suggesting  that  the  ion  might  be  able  to  tunnel  through  the  barrier  to  reach  the 
products.  We  can  make  a  crude  estimate  of  the  tunneling  effect  by  using  the  Wigner 
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expression."^  This  leads  to,,  an  enhancement  of  the  dissociation  rate  by  1.7  at  298  K.  This 
suggests  that  tunneling  at  least  needs  to  be  considered  for  predicting  the  stability  of  Ns*.  For 
lower  temperatures,  a  larger  tunneling  effect  would  be  expected. 

Lattice  Energy  Calculations 

Prediction  of  lattice  energies,  Ul  based  on  equation  (2)  requires  an  estimate  of  the 
relevant  ion  volurhes  and  gave  the  following  results: 

Ns'^Na’:  VifNs'^)  is  estimated  to  be  0.051  (+  0.020)  nm^  based  on  the  reported  ^  crystal 
structure  volume  for  N5'^Sb2Fir  by  subtraction  of  7(Sb2Fif)(table  6,  ref  30).  PCNs')  is  equal  to 
0.058  (+.0.014)  nm^  (table  5,  ref  30)  and  thus  ^(Ns'lNs')  =  0.109  (+  0.024)  nm^  leading  to  an 
estimate  for  C/l(N5'^N3')  =  142  (+  7)  kcal  mol'^ 

Ns'^Ns’:  An  upper  limit  for  V'CNs’)  can  be  estimated  from  consideration  of  y(N2)  and 
V(Ni).  The  volume  of  solid  dinitrogen,  7(N2),  which  is  dimorphous,"^^  can  be  obtained  from 
crystal  structure  data  for  the  a-cubic  and  the  hexagonal  forms  and  averages  to  V(N2)=  0.046  (+ 
0.001)  nm^  leading  to  V^fNs'^s')  =  0.155  (±  0.014)  nm^  giving  a  value  for  UlQ^s'^s)  =  129  (± 

5)  kcal  mol’’.  Alternatively,  we  can  estimate  the  volume  of  the  ions  at  the  B3LYP/6-31+G*  level 
as  described  above:  ^(Ns^)  =  0.068  nm^;  ^(Ns')  =  0.088  nm^’  ^(Ns')  =  0.065  nm^  These 
estimated  volumes  yield  C/rCNs'Tls’)  =  129  kcal  mol'*  and  C/rfTSfs'^s')  =  135  kcal  mol'*,  well 
within  the  above  error  limits.  The  corresponding  lattice  enthalpy,  AH^,  appropriate  for  use  in  the 
cycles  of  Figures  1  and  2  are  C/l(N5’*TM5')  +  2RT  and  t/L(N5‘'N3')  +  3I2RT  respectively. 

Stability  Predictions  Using  Born-Haber  Cycles 

For  a  chemical  reaction  to  proceed  spontaneously,  its  free  energy  change  must  be  zero  or 
negative.  The  decomposition  of  Ns'^Ns'  (s)  is  given  by 

Ns^Ns'  (s)  ^  2  N3  (g)  +  2  N2  (g)  (6) 

and  that  for  N5'^3'  (s)  by  reaction  (7) 

N5^3'(s)-^2N3(g)  +  N2  (g)  ^  (7) 

The  contribution  from  the  TAS  term  (entropy  change)  to  the  free  energy  for  the  reactions 
considered  in  this  paper  is  always  negative.  We  can  estimate  the  entropy  change  for  reactions  (6) 
and  (7)  as  follows.  The  entropy  for  the  solid  is  not  known  but  we  can  use  the  entropies  of 
NH4'^03'  md  NH4‘^3',  which  are  36  and  27  cal/mol  K,"*®  respectively,  to  estimate  the  value  for 
the  salts.  Given  the  well-established  entropies  of  N2  and  N3  of  45.8  and  54.1  cal/mol  K 
respectively,  we  obtain  at  298  K  for  reaction  (6)  -TAS  =  -49  kcal/mol  if  the  entropy  of 
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NEU'lMOs'  is  used  to  estimate  that  of  the  solid,  and  -TAS  =  -51  kcal/mol  if  the  entropy  of 
NHLt'^s'  is  used.  For  reaction  (7),  the  respective  values  are  -TAS  =  -35  and  -38  kcal/mol  at  298 
K.  Thus  a  significantly  positive  value  for  the  reaction  enthalpy  would  be  required  for  Ns'T'ls’  and 
Ns'^s'  to  be  stable. 

The  enthalpy  of  a  decomposition  reaction  of  an  ionic  solid  to  gaseous  decomposition 
products  can  be  estimated  from  simple  Bom-Haber  cycles  (see  Figures  2  and  3),  provided  the 
lattice  energy  (Ul)  of  the  solid,  the  adiabatic  first  ionization  potential  (IP)  of  the  anion,  and  the 
adiabatic  electron  affinity  (EA)  of  the  cation  are  known.  If  the  free  energy  change,  AGr,  of  the 
decomposition  reaction  (6)  is  positive,  i.e.,  the  decomposition  products  are  higher  in  energy  than 
the  ionic  solid,  this  value  represents  the  minimum  dissociation  energy  barrier,  but  one  should 
keep  in  mind  that  the  actual  barrier  can  well  be  higher  than  AGr.  If.  however,  the  free  energy 
change  is  negative,  spontaneous  decomposition  occurs.  The  potential  energy  curves  for  both 
scenarios  are  depicted  in  Figure  4  for  Ns'T^s'.  The  left  curve  represents  the  case  of  the  enthalpy 
of  reaction  (6)  being  negative,  resulting  in  spontaneous  decomposition  of  the  solid  to  the  first 
vibrationally  stable  intermediates,  that  is,  N3  radicals  and  N2.  The  energy  values  used  for  the  left 
curve  are  the  actual  numbers  calculated  in  the  present  study  and  demonstrate  that  there  is  no 
barrier  stabilizing  against  spontaneous  decomposition. 

The  right  curve  of  Figure  4  represents  the  scenario  of  a  positive  reaction  enthalpy  and  the 
only  condition  under  which  solid  might  be  isolated.  This  is  a  hypothetical  case,  assuming 

the  enthalpy  of  formation  of  Ns'^s'  being  about  100  kcal  mof’  more  negative  than  the  calculated 
value.  If  the  values-  for  the  IP  of  Ns'  and  the  EA  of  Ns"^  are  fixed,  the  lattice  energy  of  Ns'T^s' 
would  have  to  be  about  100  kcal/mol  higher  than  our  calculated  value  which  is  impossible  based 
on  the  molar  volumes  of  the  ions  involved.  The  relatively  small  uncertainties  in  our  AHr  values, 
-65  ±1  kcal  mol'^  for  Ns'T^s'  and  -  39  ±  9  kcal  mol'^for  Ns'TSfs',  together  with  the  large  -TAS 
terms,  allow  us  to  predict  with  confidence  that  both  compounds  cannot  be  stabilized  in  the  form 
of  ionic  solids. 

Experimental  Results  and  their  Comparison  with  the  Predictions 

All  of  our  attempts  to  prepare  solid  Ns'^'Ns'  by  the  metathesis  reaction  (8)  failed. 

Ns^SbFe'  (s)  -1-  NTNs'  (s)  — ►  MSbFe  (s)i  +  Ns^s'  (s) 


(8) 
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In  liquid  SO2,  the  Ns'^SbFe'  salt  is  highly  soluble.  Using  very  soluble  (solubility  at  -64  °C  =  3.14 
X  10'^  mol  mL’’)  CsNs,  violent  reactions  occurred  resulting  in  quantitative  formation  of  poorly 
soluble  (solubility  at  -64  “C  =  6.56  x  10'^  mol  mL'^)  CsSbFe  and  4  N2.  Using  poorly  soluble 
(solubility  at  -64  °C  =  9.23  x  10'^  mol  mL'^)  NaNs,  the  metathesis  did  not  proceed  to  any 
measurable  extent,  and  unreacted  NaNs  and  NsSbFe  were  recovered.  It  was  also  shown  that  solid 
NaNa  and  NsSbFe  can  be  mixed  as  dry  powders  at  room  temperature  and  the  Raman  spectrum  of 
the  mixture  can  be  recorded  without  evidence  of  reaction  provided  that  the  laser  power  is  kept  at 
a  low  level.  Attempts  to  carry  out  the  metathesis  in  CHF3  solution  were  also  unsuccessful 
because  of  the  low  solubility  of  the  starting  materials  in  this  solvent.  These  experiments  show 
that  in  all  cases  where  metathesis  took  place,  only  spontaneous  decomposition  was  observed  and, 
thus,  confirm  the  above  given  theoretical  predictions. 

Assessment  of  the  Performance  of  Solid  Ns'^s*  as  a  Monopropellant 

In  a  previous  publication  on  the  stability  of  Ns'^s',^  it  was  claimed  that  this  compound,  if 
it  “could  be  stabilized,  would  be  an  ideal  monopropellant  since  its  density  is  nearly  twice  that  of 
hydrazine  while  the  energy  densities  (heat  of  reaction  per  gram)  are  comparable.”  This  claim  has 
received  worldwide  attention  as  a  potential  revolution  in  rocket  propulsion.^^  Even  if  one  ignores 
the  incorrect  use  of  the  vertical  instead  of  the  adiabatic  electron  affinities  and  ionization 
potentials  in  the  previous  work,  which  as  shown  above  could  lead  to  incorrect  stability 
predictions,  the  statements  concerning  the  densities  and  energy  densities  are  incorrect.  Whereas 
the  previous  energy  density  estimate^  of  2.11  kcal/g  for  Ns’TvTs'  is  in  fair  agreement  with  our 
estimate  of  1.98  kcal/g,  the  maximum  energy  density  value  for  hydrazine  is  not  comparable  to 
this  value  at  all.  Even  if  one  assumes  its  decomposition  to  proceed  as  shown  in  equation  (9), 

3  N2H4  ^  4  NH3  +  N2  (9) 

the  maximum  energy  density  of  N2H4  would  be  only  1.17  kcal/g.^^  The  actual  value,  however,  is 
significantly  lower  because  there  is  always  some  decomposition  of  NH3  to  N2  and  H2.'‘^  Hence, 
the  energy  densities  of  and  N2H4  are  not  comparable  but  differ  by  a  factor  of  2.  The 

density  of  N2H4  is  very  close  to  1.0  g/cm^.''®  The  previous  prediction  that  Ns’lSfs'  should  have 
twice  this  density  does  not  agree  with  our  predictions.  Based  on  our  molar  volume  estimates,  the 
density  of  Ns'^Ns'  should  be  only  1.50  g/cm^.  While  the  conclusion®  that  a  hypothetical 
polynitrogen  compound,  such  as  Ns’TnIs',  would  be  an  excellent  monopropellant  for  rocket 
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propulsion  or  explosive  is  certainly  correct  and  has  been  known  for  many  years,  the  values  given 
in  reference  6  are  badly  flawed. 

Conclusions 

Our  results  are  remarkable  in  many  respects,  including  the  following  points. 

(i) :  Our  results  provide  the  first  conclusive  quantitative  predictions  concerning  the 

stability  of  Ns'T^s*  and  Ns'^s'.  They  show  that  both  compounds  cannot  be  stabilized 
in  the  solid  state.  In  spite'of  detailed  calculations,  the  previous  papers^’®’’®  dealing 
with  these  problems  did  not  use  rigorous  approaches  involving  appropriate  Bom- 
Haber  cycles,  which  prevented  any  firm  conclusions  from  being  made. 

(ii) :  Our  results  also  show  that  for  species,  such  as  and  Ns',  it  is  imperative  to  use 
adiabatic  ionization  potentials  and  electron  affinities.  Vertical  electron  affinities  and 
ionization  potentials,  obtained  by  ab  initio  calculations,  must  not  be  employed  for 
these  predictions  and  lead  to  erroneous  conclusions.  For  example,  use  of  our  lattice 
energy  values  for  Ns'^Ns'  and  Ns'lSis',  if  combined  with  the  published^  vertical 
ionization  potential  of  Ns*  (5.58  eV)  and  N3'  (2.53  eV)  and  the  vertical  electron 
affinity  of  Ns"^  (6.04  eV)  would  result  in  the  salts  Ns’^s  and  Ns'lN's'  having  large 
positive  AHr  values  with  minimum  decomposition  enthalpy  barriers  of  123  kcal  mol'^ 
and  65  kcal  mor\  respectively.  They  therefore  would  lead  to  the  incorrect 
conclusion^  that  these  salts  exhibit  extraordinary  stability. 

(iii) :  In  view  of  the  instability  of  the  Ns  radical  and  its  strongly  exothermic  spontaneous 
decomposition  to  N3  radicals  and  N2,  both,  Ns^  and  Ns ,  are  less  than  ideal  choices  for 
preparing  a  stable  ionic  nitrogefi  allotrope. 

(iv) :  Contrary  to  previous  predictions,^  Ns'T^s’  is  less  stable  than  Ns''TSr3'.  because  the 
former  contains  two  N5  species. 

(v) :  Stability  calculations  for  individual  ions  or  ion  pairs  do  not  provide  sufficient 
information  for  stability  predictions.  Rigorous  Bom-Haber  cycles,  including  realistic 
lattice  energy  estimates,  and  an  evaluation  of  the  stability  of  the  decomposition 
products  are  required  to  make  meaningful  predictions. 

(vi) ;  Our  theoretical  predictions  concerning  the  instability  of  N5''N3  and  the  use  of  the 
adiabatic  electron  affinities  and  ionization  potentials  were  confirmed  by  extensive 
experimental  studies  involving  the  low-temperature  metathesis  of  Ns'^SbFe'  with 
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different  alkali  metal  azides  in  various  solvents.  In  all  cases  where  metathesis  took 
place,  only  spontaneous  decomposition  to  Na  was  observed. 

(vii):  The  previous  conclusions®  concerning  the  energy  densities  and  densities  of 
hydrazine  and  Ns'^s*  are  not  correct. 
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Figure  Captions 

Figure  1.  Schematic  showing  the  decomposition  of  the  Ns  radical  formed  by  electron  attachment 
to  Ns*"  at  the  ZAPT2/6-311G(2d,f)  and  CCSD(T)/cc-pVTZ  levels.  The  radical  smoothly 
dissociates  to  products.  The  x  axis  corresponds  to  the  calculation  number  describing  the 
dissociation  and  the  y  axis  corresponds  to  the  energy  or  to  the  N3  -  N4  (numbered  starting  from 
the  right)  interaction  distance  representing  the  separation  between  the  N2  and  N3  molecules.  The 
large  squares  correspond  to  the  ZAPT/6-311G(2d,lf)  potential  energy  surface  and  the  circles 
correspond  to  the  CCSD(T)/cc-pVTZ  potential  energy  surface.  The  small  squares  represent  the 
N3-N4  distance  which  increases. 

Figure  2.  Bom-Haber  cycle  in  kcal  mol'^  for  the  decomposition  reaction  of  solid  Ns^s’-  The 
large  negative  value  for  AHr  demonstrates  the  unlikelihood  of  this  salt  to  exist. 

Figure  3.  Bom-Haber  cycle  in  kcal  mol"'  for  the  decomposition  reaction  of  solid  Ns'^s'.  The 
large  negative  value  for  AHr  demonstrates  this  salt  is  also  unlikely  to  exist. 

Figure  4.  Potential  energy  curves  in  kcal  mol’’  for  the  decomposition  reaction  of  solid  Ns'^s’  to 
dinitrogen.  The  left  curve  is  based  on  the  calculated  values  from  this  study  and  shows  that  a 
hypothetical  Ns’l^s’  would  spontaneously  decompose  to  2  N3  (g)  +  2  N2  (g).  The  right  curve 
shows  that  for  Ns'^Ns’  to  be  stable,  its  enthalpy  of  formation  would  have  to  be  significantly  lower 
than  those  of  the  gaseous  species  (2  N3  +  2  N2). 


NEGATIVE  B  VALUE,  KINETICALLY  POSITIVE  B  VALUE  AND  KINETIC  STABILITY 

UNSTABLE,  REPRESENTS  CALCD  VALUES  WOULD  REQUIRE  AHf  (Ng^Ns  (s))  ~  180  kcal/mol 


Activation  energy  barrier  for  decomposition 
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Table  1.  Calculated  Valence  Correlated  Geometries  Nm*.^ 


Species 

Basis 

r(NN) 

e(NNN) 

N2(X'A) 

CCSD(T)/aug-cc-pVDZ 

1.1209 

CCSD(T)/aug-cc-pVTZ 

1.1040 

CCSD(T)/aug-cc-pVQZ 

1.1005 

CCSD(T)/aug-cc-pV5Z 

1.0995b 

Expt^^ 

1.0977 

CCSD(T)/aug-cc-pVDZ 

1.1997 

u 

o 

o 

00 

CCSD(T)/aug-cc-pVTZ 

1.1840 

CCSD(T)/aug-cc-pVQZ 

1.1802 

Expt^^ 

1.18115 

CCSD(T)/aug-cc-pVDZ 

1.2072 

180.0" 

CCSD(T)/aug-cc-pVTZ 

1.1912 

CCSD(T)/aug-cc-pVQZ 

1.1876 

Expt^^ 

1.1884 

N5-(x%0 

CCSD(T)/aug-cc-pVDZ 

1.3429 

o 

oo 

b 

o 

CCSD(T)/aug-cc-pVTZ 

1.3309 

N5^(X’Ai) 

CCSD(T)/aug-cc-pVDZ 

1.1319 

165.87 

1.3413 

107.90 

CCSD(T)/aug-cc-pVTZ 

1.1159 

167.08 

1.3238 

108.48 

Expt^ 

1.10 

168.1 

1.30 

111.2 

Ns  i^Aif 

ZAPT2/6-311+G(2df) 

1.292 

121.0 

1.343 

100.1 

1.355 

109.5 

CCSD(T)/aug-cc-pVDZ 

1.422 

102.3 

1.165 

137.7 

2.576 

81.1 

Ns’-TS® 

MP2/aug-cc-pVDZ 

1.270 

125.5 

1.702 

99.4 

1.215 

107.8 

CCSD(T)/6-311+G(d)^ 

1.252 

128.9 

1.191 

1.764 

^  Bond  distances  in  A  and  bond  angles  in  degrees. 
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bCore/valence  correlation  at  the  CCSD(T)/cc-pCVQZ  level  shortens  the  bond  length  by  0.0020 
A,  resulting  in  an  estimated  CCSD(T)/aug-cc-pV5Z+CV  value  of  1 .0975  A. 

Constrained  by  symmetry  to  these  angles.  Linear  for  N3  and  N3  and  Dsh  for  N5 . 

^  The  radical  is  cyclic  with  the  symmetry  unique  atom  N1  at  the  apex.  N1  is  bonded  to  N2  and 
N5  and  N2  is  bonded  to  N3  and  N3  is  bonded  to  N4.  The  bond  distances  are  in  the  order;  r(Nl- 
N2)  and  r(Nl-N5);  r(N2-N3  and  r(N4-N5);  and  r(N3-N4).  The  bond  angles  are  m  the  order 
<N2N1N5;  <N1N2N3  and  <N1N5N4;  and  <N2N3N4  and  <N5N4N3. 

*  The  transition  state  for  the  anion  decomposition  is  cyclic  with  the  symmetry  unique  atom  N1  at 
the  apex.  N1  is  bonded  to  N2  and  N5  and  N2  is  bonded  to  N3  and  N3  is  bonded  to  N4^'lje 
bond  distances  are  in  the  order  r(Nl-N2)  and  r(Nl-N5);  r(N2-N3  and  r(N4-N5);  and  r^3-N4). 
The  bond  angles  are  in  the  order  <N2N1N5;  <N1N2N3  and  <N1N5N4;  and  <N2N3N4  and 
<N5N4N3. 
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Table  2.  Calculated  Valence  Correlated  Harmonic  Vibrational  Frequencies  (cm'l). 


Mode 

Mode 

Mode 

Method 

Number 

Symmetry 

Description 

CD 

Expt. 

Na 

CCSD(T)/ 

1 

Na  stretch 

2319 

233  f 

aug-cc-pVDZ 

(2358.6) 

CCSD(T)/ 

aug-cc-pVTZ 

1 

Na  stretch 

2340 

CCSD(T)/ 

1 

c^g 

Na  stretch 

2354 

aug-cc-pVQZ 

CCSD(T)/ 

1 

c^g 

Na  stretch 

2359” 

aug-cc-pV5Z 

Ns 

MP2/cc-pVTZ 

1 

71 

bend 

594 

457® 

2 

Qcr 

c> 

symmetric 

stretch 

1434 

1320® 

3 

Ct« 

asymmetric 

stretch 

1620 

1645® 

Ns' 

CCSD(T)/ 

aug-cc-pVDZ 

1 

TC 

bend 

578 

626*^ 

2 

<^g 

symmetric 

stretch 

1281 

1335® 

1986^ 

3  . 

CTu 

asymmetric 

stretch 

1999 

Ns' 

CCSD(T)/aug 

-cc-pVDZ 

1 

ei’ 

NN  stretch 

1202 

2 

ai’ 

NN  stretch 

1141 

3 

ea’ 

NN  Stretch 

1078 

4 

bend 

1001 

5 

e2” 

bend 

739 

MP2/ 

1 

ei’ 

NN  stretch 

1143 

aug-cc-pVDZ 

2 

ea’ 

NN  stretch 

1123 

3 

ai’ 

NN  Stretch 

1094 

4 

ea’ 

bend 

1039 

5 

ea”  ' 

bend 

740 

N5^ 

MP2/ 

1 

ba 

NN  str 

2260 

2267® 

aug-cc-pVDZ 

2 

ai 

NN  str 

2179 

2206® 

3 

ba 

NN  str 

1199 

1064® 

4 

ai 

NN  str 

881 

870® 

5 

ai 

bend 

669 

671® 

MP2/ 

aug-cc-pVDZ 


ZAPT2/6- 

311+G(2df) 


6 

b2 

bend 

431 

4148 

7 

a2 

out-of-plane 

bend 

.426 

474® 

8 

bi 

bend 

374 

421® 

9 

ai 

bend 

Ns'-TS 

179 

204® 

1 

b2 

NN  stretch 

1649 

2 

ai 

NN  stretch 

1572 

3 

ai 

NN  stretch 

1244 

4 

ai 

NNstretch 

905 

5 

b2 

bend 

848 

6 

b, 

bend 

684 

7 

a2 

out-of-plane 

bend 

591 

8 

b2 

bend 

182 

9 

ai 

NN  stretch 

Ns  (cyclic)  ^Ai 

862i 

1 

b2 

NN  stretch  ■ 

1221 

2 

ai 

NN  stretch 

1147 

3 

ai 

NN  stretch 

1064 

4 

b2 

NN  stretch 

1033 

5 

ai 

NN  stretch 

1002 

6 

ai 

bend 

843 

7 

b2 

bend 

798 

8 

a2 

out-of-plane 

bend 

729 

9 

bi 

bend 

363 

*  Ref.  33.  Value  in  parentheses  is  the  harmonic  value. 

”  Core/valence  correlation  at  the  GCSD(T)/cc-pCVQZ  level  of  theory  increases  the  harmonic 
frequency  by  9.2  cm'\  resulting  in  an  estimated  CCSD(T)/aug-cc-pV5Z+CV  value  of  2368 

cm  . 

"Ref.  38 

^  Lamoureux,  Dows,  D.A.;  Spectrochim.  Acta,  Part  A,  1975,  31, 1945.  From  N3  in  a  KI  lattice. 
"  References  40  and  41. 

^Reference  35. 


^  Reference  2. 
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Estimated  frozen  core,  complete  basis  set  energy  obtained  from  Eq.(l)  using  the  CCSD(T)/aug-cc-pVxZ  (x  =  D,  T,  Q)  energies.  For  the 
sake  of  comparison,  the  CBS  dissociation  energy  of  N2  with  the  aug-cc-pVTZ  through  aug-cc-pV5Z  basis  sets  is  226.63  kca/mol. 


